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We study the dissipative thermodynamics of vortex lines in layered superconductors within a simple 
string model in the dilute limit of negligible vortex interactions and compute the specific heat C v 
in presence of arbitrary dissipation. The interplay of dissipation, inertia and elasticity is shown to 
control the qualitative thermodynamical behavior and their relative amount determines two very 
distinct regimes for the specific heat. In the dissipation dominated case we find a behavior C v ~ yT 
for a large interval of temperature below T c . 
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I. INTRODUCTION 
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As well known, type-II superconductors in a magnetic 
field exhibit an intermediate mixed phase characterized 
by vortex condensation and partial magnetic flux pene- 
tration through them. The dynamics of the latter cru- 
cially drives the electromagnetic behavior, and has been 
extensively studied in the last yearsfcl. It has been shown 
that under appropriate assumptions, in particular for 
small magnetic field and vortex density, the elementary 
vortex excitations in thin layered samples with trans- 
verse magnetic field are flux lines crossing each layer only 
once, joining pancake vortices forming in each layer and 
moving in the xy.-plane. A well established and studied 
continuum modellil exist for these string-like excitations 
and it has-,been succesfully applied 4,o analyze both the 
dynamicaBa and therodynamicalEiQ behavior. 

In the following we study the specific heat of vortex 
lines within a simple model which can be treated exactly 
and captures all the essential features of vortex dynam- 
ics relevant in layered superconductors. Dispersive and 
magnetic effects will be neglected, whereas dissipation 
(rj), inertia (/x) and elasticity (e) will be taken into ac- 
count exactly in the general case of a finite number N 
of layers. We shall call d the inter-layer distance, N the 
number of layers, L the thickness of the film and A its 
area. The continuum limit N — > +oo will be analyzed as 
a particular limit in which the usual model holds. The 
problem has already been addressed in the latter context 
by Blatter and IvlcvlLfor small magnetic fields and by 
Bulaevskii and Maleyu for large magnetic fields. 

One of the central issues of this paper is to present a 
peculiar and unambiguous relation between the value of 
the ratio / = and the thermodynamical properties 
of the vortex line. We show that there is a tempera- 
ture T a — f j^-^r, typically of the order of a few to A, 
which separates the universal /-independent low temper- 
ature already discussed by Blatter and IvlevO from an 
/-dependent higher temperature (but still below the crit- 
ical temperature) behavior. In the relevant overdamped 



case in which / <C 1, we find that in this regime the 
specific heat is proportional to the square root of the 
temperature. This should hold up to the critical temper- 
ature, provided that the model is still appropriate. In 
Section IV we summarize our results; we compare them 
with some previous results appeared in the litterature 
and discuss their physical relevance. In section V we dis- 
cuss the effect of including the Magnus force. We find 
essentially the same results as before, in particular the 
same behavior with temperature, and even the numeri- 
cal values of the coefficient are not very different. 

The fundamental characteristics of the vortex lines are 
quite simple in the approximation we are considering; in 
particular, the self-energy stemming from the interaction 
between neighbor pancake vortices essentially reduces to 
a local elastic term, and dissipation seems to be reason- 
ably well described by an ohmic viscous term. Both the 
elastic modulus e and the friction coefficient r\ can be the- 
oretically computed, but the mass density /x is still object 
of discussion and can be only hardly estimated because 
of the wide class of contributions it gets depending on the 
particular setting, and its determination is still a central 
issue of vortex physics, even if there is strong belief for 
it to be negligibly small with respect to dissipation^. 

Since each point of the string has two-dimensional dy- 
namics, we parameterize the problem in cylindrical co- 
ordinates, calling q(i, z) the xy position of the string at 
given z. The equation of motion of the vortex line is 



/iq + nq - eq 



. 



(1) 



Because of the finite size of the system in the z direc- 
tion, it is necessary to impose the additional Neumann 
boundary conditions 



q'(t,0)=q'(t,Z)=0 



(2) 



In this way, no energy-momentum flow is allowed across 
the end-points of the string. 

The quantum-mechanical version of this dissipative dy- 
namics is reached by means of Caldeira and Leggett's 
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formalism!, in which dissipation is introduced by cou- 
pling the system to a thermal bath of harmonic oscil- 
lators with a continuous frequency spectrum satisfying 
the requirement to reproduce a simple velocity depen- 
dent ohmic viscous term in the effective equation of mo- 
tion. The Euclidean effective action of the open system, 
relevant for thermodynamics, is obtained by integrating 
out the degrees of freedom of the bath, and reads in the 
continuum limit 
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S2„ = sin Si 

7T \2NJ 
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-fiq (r,z) + -eq (r, z) 
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q(T, z)-q(T\z) 
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It is convenient to extend q(t, z) symmetrically to neg- 
ative z , and expand the obtained even function in Fourier 
modes of fundamental frequency v = ^ along z by means 

of q(i, z) = q (t) + V2~Y,n=o cos vnz; in this way, 

the boundary conditions arc automatically satisfied. The 
Euclidean action then splits into a sum of harmonic os- 
cillator contributions encoding the vibrational modes of 
the string plus a zero mode representing the translational 
motion of the center of mass. 

In order to handle the thcrmodynamical Euclidean 
path-integral, we introduce the Matsubara frequency 
ui = -0- and expand the periodic configurations q n (r) 
involved in the path-integral for the partition function 
in Fourier modes with fundamental frequency u> along 
imaginary time, taking q„(r) = -±= Y.t=-oo Qnke 1 ^ ■ 

The dissipative term, once Fourier-transformed, will con- 
tribute a peculiar \uk\ term, remnant of the derivative 
nature of ohmic dissipation. 

In Fourier space, the effective action is 



N-l +oo 
n— k— — oo 



■A\w k \] |q„ fc | 2 . (4) 



M = uL is the total mass, A = r\L the total friction 

coefficient and SI = . the characteristic vibration 
V ^ 

frequency; we have used the notation u>k = tok and Sl„ = 
fin. Since there is only a finite number N of layers, we 
have cut-off the possible modes at that value. 

A more natural way back to the finite N case is to 
write a discrete action in terms of the single pancake 
vortices coordinates x;(£) = q(t,ld). One then obtains 
the action for N pancake vortices of mass fid and friction 
coefficient 7/g?, each harmonically coupled to its nearest 
neighbors with elastic modulus %■ The decoupling of 
these degrees of freedom can be achieved performing a 
change of variable analogous to the Fourier expansion of 
the continuum case, defining q„(t) from x;(i) through 
= q a (t) + \^J2 n =i Jin(t) cos (vn(l + ±)d) . One 
then obtains the action (|4|), but with a modified fre- 
quency spectrum given by 



(5) 



Obviously, in the limit of N — ► +oo, the lattice version 
of the theory reduces to the continuum formulation, with 
the same frequency spectrum Sl„ = Sin. 

In the following, for the sake of clearness we will con- 
centrate on the basic case of absence of pinning of vortex 
lines by impurities and defects, and focus on the effects 
of elasticity and dissipation (the effect of the Magnus 
force will instead be considered in Section V). However, 
the treatment of this important generalization could be 
easily faced in exactly the same way; introducing a sin- 
gle columnar pinning center through a harmonic con- 
fining potential with elastic modulus k and characteris- 
tic frequency Q p = yfjj-, the only novel feature would 
be the modification of the spectrum of the modes to 

n: = 
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ii. PARTITION FUNCTION 

We have seen that the dynamics of the vortex line is 
encoded in a zero mode describing its translational mo- 
tion and N — 1 harmonic modes with frequency spectrum 
Q n describing its vibration. Consequently, the partition 
function, 



Z(f3) = Tre 



-0H 



Dq(r, z)e 



-S(q(T,*)) 



(G) 



q(0,2)=q(/3,. 



will factorize into the product of the partition functions 
of the modes. Thus, all that we need to know for our 
aim is the dissipative thermodynamics of the zero-mode 
degree of freedom and the remaining harmonic modes. 
The path-integral for the latter can be evaluated quite 
easily in both cases since it is Gaussian, and reduces to 
a product of ordinary integrals in Fourier space. 



A. Zero-mode 



For the zero mode, we get 



Z{j3,h)=N{p) J] // dq ok dq: k e-^( M ^ +A ^h^ 
&=()•'•' 

+oo / . \ -2 

=M'mu (i+i- h ) . (7) 



k=l 



where 7 = -^j = 2 is the characteristic frequency related 
to dissipation and A is the area. The normalization factor 
is fixed by requiring that in the limit A — > the partition 
function reduces to the free one Z {j3) 



tt^A. It follows 
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(8) 



Unfortunately, the infinite product representing the ef- 
fect of dissipation is divergent. We then have to introduce 
a frequency cut-off uj c above which the effect of dissipa- 
tion is assumed to be small, or at least not adequately 
described by a simple ohmic viscous term in the equation 
of motion, and truncate the product at k c (ui) = —. Us- 
ing t he i nfinite product representation of the T-function, 
Eq. (A2), we obtain for k c (ui) — > +00, 



k c (u) 

n 
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(9) 



it cap_be obtained using the Drude model 
3 nB|n|, , 



The same result 

as regularizationElEil, where the frequency spectrum of 
the thermal bath is cut-off at lo c dividing it by a factor 
1 + (^-) 2 ; the cut-off frequency can then be related to a 

microscopic relaxation time r c by 10 c = — . 

The final form for the partition function is conveniently 
written in terms of the function 

f 0) = -Lz5-z e *r(z) , (10) 



/27T 

which goes to one for la rge arguments away from the real 
negative axis (see Eq. ( |A5|) ); one gets 



Z(f3, 1 ) = A^f 2 

Z7T 

with a zero-point energy 
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£,(7) = - 1 
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(11) 



(12) 



Quite interestingly, we learn that, since the partition 
function depends only on — apart from an irrelevant mul- 
tiplicative constant and the zero-point energy, the effect 
of dissipation on the thermodynamics of the zero-mode 
is only to modify the temperature scale. This could have 
been expected from dimensional analysis, since the only 
energy scale we can construct from M and A is 7 (we do 
not consider lo c since we expect dissipation to manifest 
in a universal ohmic way in the limit of small relaxation 
time r c , situation in which the only role of w c should be 
to renormalize the energy of the open system) , and the 
only dimcnsionlcss quantity that can enter the partition 
function is the reduced temperature — . 

Notice that in the limit of strong damping ^ ^> 1, that 
is equivalent to the small mass limit fj,— > since — = 
^p, we are left, apart from the zero-point energy, with a 
trivial constant par tition function Z[((3,"f) independent 
of //; using Eq. rfA9) 



£>1 27T 
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B. Harmonic modes 

For the other modes, we proceed exactly in the same 
way and the partition function is given by 



+00 
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The requirement this to reduce to the pure harmonic os- 
cillator result Z (j3,tt n ) = (2sinh^) -2 in the limit 
A — > fixes, using Eq. (Al), the normalization coeffi- 
cient to 

= J 2 ■ (15) 

As before, the dissipative term is responsible for a diver- 
gence in the correction factor accounting for the damping 
of the system; introducing the same frequency cut-off, the 
product can be carried out by factorization, obtaining 



n 

fc=l L 



1+-T- 



u> k 



1 

k 2 



(^) 2 r(^)r(i^) 



(16) 



where £„ = y f2 2 — ^- . Again, the same result can be 

obtained using the Drude model as regularization00, 
allowing a microscopic interpretation of the cut-off. 

The partition function can then be written, using again 
the T-function defined above, as 



z(j3,n n n) = f 2 



2 i£n 



-/3£ (n„, 7 ) 



(17) 



with a zero-point energy given by 



f (n n)7 ) = l(l + ln^)-i^ln| , 

7T V il n J 7T i- It,, 



1 



(18) 



In the strongly dissipated regime, which means ^ n, 
and successively in the limit of low temperature with re- 
spect to damping, corresponding to ^ 3> 1, th e part itio n 
function is seen, using respectively Eqs. (A4) and (A5), 
to reduce to 



I27 

n 2 

n 
'JUJ 



f 2 ( ^l] e -P£o(~f) 
\7<<V 



-0£o(rr) 



(19) 
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The last extreme limit can be seen as a small mass 
limit fi — > 0. Since the variable = is in- 

dependent of fi, the resulting behavior Zi((3, f2 n , 7) no 
longer depends on fj, apart from the zero-point contribu- 
tion £0(7), which is found to be the same as the one for 
the zero- mode, Eq. (|J). A gain, this is obvious from 
dimensional analysis. The only dimensionless quantity 
independent of fi we can form from the three variables 
w, 7 and f2 is indeed entering the resulting behavior; 
instead, the /u-dependent variable ^ = ^ disappears in 
the limit. Having eliminated one of the scales, we fall 
in a regime in which, as for the zero-mode, dissipation 
manifests itself only in the temperature scale. 



C. Vortex-string 

In order to compute the partition function of the 
string, all that we have to do is multiply those relative 
to the modes that we have just computed, that is 

N-l 



n=l 



"--2 / 2 



p2 I 2 



-/3i5 (n, 7 ) 



(20) 



The zero-point energy is the sum of those of the modes 



E f in (n,~f) 
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(24) 



Notice also that 



The last sum is finite also for N — > - 
it vanishes in the free case 7^0. 

In the limit in which each of the N modes is in the 
strong damping regime and at low temperatures with re- 
spect to friction, meaning ^ ^> N and ^ 3> 1, and cor- 
responding to a limit of vanishing mass /1 — > 0, the parti- 
tion function first factorizes essentially into the product 
of two functions depending respectively on the combina- 
tions & = if and $ = (f ) 2 and then simplifies 
to a function Zi vor {(3, f2, 7) which, apart from the zero- 
point energy, depends only on the /i-independent variable 
In fact, 



Z vor {[3,Cl,~f) 



— ► A^f 2JV 
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7W 
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(25) 



JV-l 



S (fi l7 )-f (7)+ ^5 (fi„, 7 ) 



(21) 



The first observation we can do is that for fixed w, 
7 and f2, the factor Z((3, £l n ,l) in the product over the 
modes in the partition function goes to one for n — ► +00, 
and the latter is thus a well defined function of these vari- 
able for any N, since the product can be shown to con- 
verge. The meaning of this observation is that the vari- 
ables we have chosen represent the true physical scales of 
the problem, and increasing TV just add higher frequency 
modes that are more and more suppressed on the fixed 
scale we consider. 

Conversely, the zero-point energy increases indefinitely 
with N. Using the lattice regularization and the results 
derived in appendix B, the total zero-point energy can 
be recast in the following form 



E (Sl, 7) = E odiv (Q, 7) + E ofin (Q, 7) 



(22) 



with 

E odlv (n,j) 



N [AN \ n 1 7 2 
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where now the zero-point energy has simplified to 



£ ( 7 ) = N£ (>y) = 

TT 



1 



111 

7 



(26) 



Again, we see that the product in the partition func- 
tion is well defined for any N since the n-th factor goes 
quickly to one for n — ► +00; the zero-point energy grows 
indefinitely if one add modes increasing N, as in the gen- 
eral case. 

In the following, we shall distinguish between two 
regimes in the thermodynamics of the vortex line: the un- 
derdamped one, for which % <C N, and the overdamped 
one, for which ^ ^> N. 



III. SPECIFIC HEAT 

Having computed the partition function of the vortex 
line, we can easily compute its specific heat. The latter 
can be defined in presence of dissipation as the derivative 
of the mean energy of the open system with respect to 
temperature, considering dissipation just as an additional 
mechanism for exchanging energy with the thermal bath 
and modifying the response of the system under heating. 

In fact, within the Caldeira-Leggett formulation, the 
partition function of the complete system factorizes into 
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the partition function of the bath and that of the open 
system with an effective action containing the influence of 
dissipation. Technically, this is implemented with an ap- 
propriate shift in the bath's oscillators coordinates which 
are just dummy variables in the path-integraEHi3. Thus, 
in this scheme, the total heat capacity of the system is 
the sum of the heat capacity of the bath alone plus the 
one of the vortex lines. 

One of the interesting features of the specific heat is 
that it is independent of the rather unknown parameter 
uj c entering the zero-point energy. Also, from our previ- 
ous discussion, we expect it to remain well behaved in 
the ji — > limit. 



N-l 



-> 27V 
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(30) 



In the final extreme limit, the specific heat reduces to 
a function C;„ or (/3, fi, 7) which depends only on the fi- 
independent variable S?. 



A. General expression 

The specific heat of a single vortex line is 
& 2 



dTd/3 



lnZ str (/3,n,7) 



(27) 



From now on, we will work with the reduced specific 
heat, dropping a factor kg. The usual one is obtained 
by multiplying the latter by fee and by the density of 
vortices in the sample. 

The result for the specific heat is easily computed as a 
sum over the contributions of the modes and is found to 
be 

N-l 

C vor ((3, 0, 7) = C{fi, 7) + E n »> 7) 



n=l 



2 _ 




where 



*'(*) = ^3 lnf(*) 



(28) 



(29) 



and goes to zero for lar ge arguments away from the real 
negative axis, Eq. (A6). At high temperature, the spe- 
cific heat saturates at the value 2N — 1; all the TV — 1 
harmonic modes contributes 2, whereas the zero- mode 



only 1, as seen from Eq. (A8). 

In the limit of strong damping for each of the N modes, 



% ^> A, we have 



— 2. At low temperatures 



with respect to friction, % ^> 1, the specific heat simpli- 
fies to 



B. Behavior in various regimes 

As shown, the variables which naturally enter the ther- 
modynamics of the string have the dimensions of an en- 
ergy and are the Matsubara frequency ui = 27rfcsT which 
is proportional to the thermal energy, the characteristic 
frequency Q — °^ ^ nc string and 7 = ^ which is 

the characteristic energy of dissipative processes. These 
variables represent the dependence on temperature and 
damping on a physical energy scale that depends on 
L = Nd and thus on N, the number of layers. 

It is convenient to define the reduced temperature t 
and the dimensionless damping parameter a according 
to 



/ = - = l\'-dk R T =— — 



e 

7 N rjd 

n = 



N T n 



efj, 



N 
^1 



(31) 
(32) 



The /i-independent variable arizing in the /i — > limit is 
then simply 



at 



u>7 

ft 2 " 



2N 2 r)d 2 



T 



(33) 



The meaning of a possible continuum limit A — > +00 
can now be better elucidated. In order to compare situ- 
ations with different TVs, we work with the A-dependent 
variables a and t, so that adding modes does not change 
the contributions of those that were already present but 
constitute only a correction. The continuum limit will 
then be relevant as a reliable approximation to the finite 
A case whenever the modes with n > N give a negligible 
contribution to thermodynamics. As we will see, this will 
be true for low temperatures. 

As we already pointed out, the thermodynamics of the 
vortex line clearly exhibits two distinct and very different 
regimes. We shall call the underdamped regime the case 
in which a <C A, and the overdamped regime the case in 
which a A, and analyze them separately. 

The global behavior of the specific heat is shown in Fig. 
0. There are two regimes according to the value of a. For 
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a <C N , the underdamped case, the specific heat depends 
only weakly on it and the effect of increasing TV is to raise 
the temperature at which the saturation to the asymp- 
totic value 2N — 1 begins. For a ^> N, the overdamped 
case, the specific heat acquires a strong dependence on 
a. Increasing the temperature, the specific heat raises 
very quickly to the value N — 1, going approximately 
like a square root, and continues to grow very slowly and 
almost linearly towards its high temperature saturation 
value 27V - 1. 

In the limit of very small mass, the specific heat stop 
its raising at N — 1 instead of 2N — 1. This is due to the 
fact that the extreme damping has killed almost all the 
kinetic part of dynamics, corresponding to an asymptotic 
specific heat equal to instead of 1 for the center of mass 
particle degree of freedom (which gets completely killed) 
and equal to 1 instead of 2 for each of the harmonic modes 
(for which the kinetic and potential part of the energy are 
equal in the free case). Moreover, if damping is strong 
enough with respect to inertia, the specific heat depends 
only on the /^-independent variable at. 




FIG. 1. The specific heat C vor (t,a) of the vortex-string as 
a function of the reduced temperature t. The solid lines refers 
to a = 0, 3, N, 2N, A.5N, 12N. For a = A.5N, 12N, the lim- 
iting strong damping scaling behavior with (dot-dashed lines) 
and without (dashed lines) the linear correction (term y ^ in 



C vor (t, a) 



(1/1 

sU + r 1M<1 



-t, l«f«JV 

o 



(34) 



2N — 1 , t > N 



The behavior for t -C 1 directly follows from the asymp- 
totic behavior of the w- function, Eq. (Ah), whereas the 
one for 1 -C t <C -V is obtained approximating the sum 
in the specific heat with an integral. 

At low temperatures, t <C N, the underdamped case 
further splits into three sub-regimes. For values of a be- 
low a first critical value a c ± ~ 0.4, the specific heat pro- 
gressively pass from the free shape to a linearly growing 
one; moreover, any infinitesimal amount of dissipation, 
a > 0, causes the specific heat to start from zero in- 
stead of one, because of the center of mass contribution. 
For a beyond a second critical value a c i — 1.5, the spe- 
cific heat starts to deviate from the linear shape, getting 
more and more convex at low temperature. Finally, for 
a between the two critical values, there is no substan- 
tial dependence of the specific heat on friction, and the 
response to heating is approximately linear in the tem- 
perature in the whole continuum region t -C N. 

The two critical values of a can be obtained by re- 
quiring the matching of the behaviors for t <C 1 and 
1 <C t <C N. They satisfy thus the quadratic equation 



which yields 



a c i,2 
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(35) 



3T V3 



(36) 



Eq. 



is seen to be an increasingly good approximation. 



in agreement with the values extracted from the plot. 



2. Overdamped regime: a> N (I <C 1) 



1. Underdamped regime: a N (I ^> 1) 

In this regime, we use Eq. (pffi), which depend on both 
a and t. Notice first that in the high temperature limit, 
t ^S> N, the specific heat tends to its maxmimal value 
2iV — 1 and is almost constant. Next consider not too 
high temperatures, t -C N. The contribution of the 
modes with n > N in the sum giving the total specific 
heat would give a negligible contribution (since the argu- 
ments of both the ^-functions would be large) , and the 
continuum limit is a good approximation. We can thus 
take the limit N — > +oo keeping a and t fixed. In this 
way, the asymptotic behaviors of the specific heat in the 
extremes of this range of temperature can be computed 
and finally 



In this regime, we can use Eq. (|3ti|), which depend only 
on at — ^r. At high temperature, at ^> N 2 , the specific 

heat grows linearly from the value N — 1 with a slope ^ 
until at ~ a 2 , where it starts saturating to its asymptotic 
value 2A^ — 1. For low temperature limit, at <C N 2 , the 
contribution of the modes with n > N in the sum giving 
the total specific heat are again negligible, and we can 
take the continuum limit N — > +oo keeping at fixed. In 
this case, the asymptotic behaviors of the specific heat at 
the extremes of the low temperature region can be found 
and finally the behavior is the following 
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C vor (t, a) 
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2C Vat , 1 « at « JV 2 

N , , TV 2 

iV - 1 + — * , TV 2 < at < a 2 = 



2iV - 1 . at > a 2 



TV 2 

7T 2 / 2 



(37) 



where C = 0.490 is the constant given by the integral 
(All) quoted in appendix A. As before, the behavior for 
at <C 1 directly stems from the asymptotic behavior of 
the ^'-function, Eq. (At), whereas the one for 1 <C at <§; 
iV 2 is computed approximating the sum over the modes 
with an integral. 



(40) 

For reasonable AT, ranging from 10 2 to 10 4 and corre- 
sponding to a thickness L of the sample between 0.1 pm 
and 10 /urn, T s can go from 10 misT to 1 K, indeed repre- 
senting the accessible temperature scale for the problem. 
In the following, we will assume T <C j^jt^ji m order 
to satisfy the condition 3- ^> 1 which has been used in 
deriving the behavior in the overdamped case. Using the 
estimate given at the end of this section for /z, this means 
temperatures below 10 5 K, and does not constitute any 
restriction. 

In terms of T s , the reduced temperature is given by 



t ■■ 



e/Z T 



r]d T s 



(41) 



The reduced temperatures that are important for the 
shape transitions in the specific heat are recognized to be 
t = 1, N in the underdamped case and t = — , a in the 
overdamped one. The important question that we shall 
now address is whether some of these can be relevant at 
the superconductivity temperature scale. 



IV. DISCUSSION 

The available experimental and theoretical work on 
superconductors and vortex dynamics does not allow 
for a precise knowledge of all the parameters entering 
the description of the problem. In particular, whereas 
the damping and elasticity coefficients can be computed, 
the mass density is ambiguous. We will thus focus on 
J = as the fundamental unknown quantity carrying 
the dependence on ji. 

The value of rj can be estimated atj-Lpw temperaure 
using the Bardeen-Stephen expressiontj, whereas e is 
known from the microscopic theory; in Gaussian units 
(and recovering h), 



$ 2 



\47rA / k£ 



(38) 



where A is London's penetration depth, £ the xy coher- 
ence length, k = y/jfi tne anisotropy ratio, the nor- 
mal state resistivity, and <5 = || the flux quantum. 
For YBCO films, we can take the typical values d ~ 

12 A, A ~ 1400 A, £ - 15 A, k ~ | and p N ~ 100 pVL-cm, 
yielding 



T]~ 3.0 10" 



_Erg ■ s 



e~ 3.4 10 



-7 Er 9 



(39) 



Observe that since T s is known but /x is not, the tem- 
perature scale entering t is ambiguous. We are now able 
to look closer to the order of magnitude of the transition 
temperatures in the two regimes we have studied, tak- 
ing as reasonable temperature scale T s (T = is thus 
small). 

In the underdamped case, a = ^"^f^ <C N, the im- 
portant reduced temperatures where found to be t& = 1 
and ts = N, corresponding to the temperatures Ta — 
§T S < T s and T B = aT s ~ T s . Thus, we conclude that 
the relevant regimes for the specific heat are in this case 
the first and second of Eq. (|34|) , that is 



N ( 1 tt 2 \ T T s 




N 



(42) 



In the overdamped case, a = -^^L ^> jV, we had 



instead tj± — ~t ts = ~~ and tc — a, corresponding to 
the temperatures T A = ±T S < T s , T B = NT S > T s and 

2 

Tc — jj-T s ^> T s . Thus, the relevant regimes for this 
case are the first and second of Eq. (|37|) , that is 



Ntt 2 T 



T 

T« — 

N 



(43) 



Multiplying the results ( |42| ) and (43) by ks and the 
density of vortices in order to obtain the true specific 
heat per unit volume, and expliciting all the variables, 
our final result can be written as 



In order to get contact with realizable temperatures 
in the framework of superconductivity, we define the fi- 
independent temperature 
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(44) 
/< 1 
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(45) 



2 fcs r/L 2 ' 77c? 

In a previous work by Blatter and IvlevB the low tem- 
perature regime for small magnetic fields was considered 
using the usual continuum description; this corresponds 
to the limit N — > 00 in our description and temperatures 
T -C T a . In the anisotropic and elastic limit, Eq. (33) of 
that work for the specific heat reduces to 



C,, 



2 B fc| 77 
3$ ft e j. 



1.2 



(46) 



Our finite N case can be recovered with the identifica- 
tions k z — > ^n, dfcj — > (iv = ?•), meaning k zm i n = ~ 
and 



n=l 



G 



(47) 



In this way, one matches the linear and /^-independent 
results of the first and third row of (Q) , which are rele- 
vant for T -C T Q and coincide in the continuum limit on 
wich we are focusing. 

Note that out result Eq. (^) , case T ^> , would be 
modified by the introduction of a lower cut-off k zm i n in 
the following way 



sf 



dzz 4 &(z 2 



N- 



We took k. 



but the behavior C vc 



would hold also for other possible k z 



(48) 



provided 



A- 



< 



The behavior for large magnetic fields has been stud- 
ied by Bulaevskii and MaleyQ. In this case, the specific 
heat is still linear in the temperature, but has a different 
dependence on the magnetic field and the microscopic 
parameters of the vortices, and the result can not be 
expressed as a function of only u, r\ and e; this signals 
that in this regime the vortices do not retain their line 
structured and other kind of configuration become im- 
portant. 



Summarizing, for T <C T Q , and in the continuum limit 
N — > 00, we fincLthe universal behavior already studied 
in the litteratureEl. For T T Q and in the whole range of 
superconductivity temperatures we find instead a differ- 
ent behavior depending on the value of the parameter /. 
In the underdamped case / 3> 1 the specific heat is lin- 
ear in the temperature and depends only on inertia and 
not on friction, whereas in the overdamped case I <C 1 it 
goes like the square root of the temperature and depends 
only on friction and not on inertia. 

It is important to observe that for N ranging from 
10 2 to 10 4 , the temperature T Q is very small and varies 
from 1 ixK to 10 mK. This suggest that an interesting 
experimentally observable regime should be T > T Q . 

As an example of theoretical estimate, one can approx- 
imate the mass density fi with its electronic and elecf 
magnetic contributions. These are found found to be 



fJ-el 



-m e kF , /Ite 
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(49) 



For YBCO films, the Fermi momentum is kp = 

o 

0.5 A^ 1 , and one obtains u e i — 2.9 10~ 21 (i em ~ 
12 10-22 <£__ This wou id mcan J ~ o.l, which does not 
corresponds clearly to any of the two damping regimes. 
Thus, if one uses this theoretical estimate, the ther- 
modynamics should lie somewhere inbetween the under- 
damped and overdamped regime that we have consid- 
ered in mope detail. However, trusting other conventional 
arguments^, the overdamped case should be the relevant 
one. An experimental measure of the specific heat for 
T> T„ could thus provide specific informations on this 
important issue. 



V. EFFECT OF THE MAGNUS FORCE 

In this section we will discuss the influence of a possible 
Magnus force on the thermodynamics and show how our 
results for the specific heat generalize to the case in which 
both the Magnus effect and friction are important. 

The Magnus force corresponds to a term 



<5*q 



(50) 



in the equation of motion (Q) of the vortex line, where 
*q = q x z is the dual of the xy position. This term can 
be accounted for in the thermodynamics by adding in the 
Euclidean effective action (|J) the term 



dr—Sq *q 
2 



(51) 



dz / 
'0 Jo 

Actually, the Magnus force is rather controvertialEfl. In 
the low magnetic field limit, by estimating S from the 
number density of the superconducting fluid, one would 
get a value of the same order of 77 as in eq. (B9) . 
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Accordingly, the partition function for each of the 
modes of the vortex line engoes the following modifica- 
tion 
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(52) 



We have introduced the complex damping frequency 
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(53) 



This leads to the following modifications in the expres- 
sions for the partition functions 



r 



rillf lr 



■ 1 + if 
p2 / 2 ^ 



(54) 



The constant appearing in the overdamped case 
now depends on the phase <fi and is again given by an 
integral 



-2i0 



dzz 4 *'( e -^z 2 ; 



(59) 



Since the relevant asymptotic behavior of the integrand 
for z — > is independent of 0, we expect C$ to depend 
only weakly on it. In fact, the to extreme values corre- 
sponding to r\ ^ 0, S — ((f) — 0) and 77 = <5 = ^ 
(^ = 5) are given by the integrals ( All ) and ( A12 ) of 
appendix A, 



C = 0.490 , Cjl = 0.346 



(60) 



For arbitrary 0, is of the same order of magnitude 
and can be computed numerically. 

Observe finally that in the overdamped case, in the 
sense explained above (with respect to ?y e ), the specific 
heat depends on the square root of the temperature and 
on the effective damping coefficient r\ e = y^r/ 2 + 8 2 . In 
the ultra clean limit ir\ e = S) this was already noticed 
bye Blatter and Ivlevu, who also found C v ~ VT in this 
case. 
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where 



e« = - y , t = yw n -^. (55) 

All the analysis done for the dissipative case goes 
through without any major difficulty. Underdamped and 
overdamped cases now means with respect to the mod- 
ulus rj e of the complex friction coefficient i) — r\ e e % ^ for 
which 



Ve = v V 2 + $ 2 i 4> — arctan — 

77 



(56) 



The result ( p4[ ) generalizes to 

( ir B kg ( r]L 6/J77 
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APPENDIX A: PROPERTIES OF T-FUNCTIONS 

In this appendix, jue recall some formulae involving T- 
functions and otheilij. First of all, there are the infinite 
product representations of the hyperbolic sine and the 
T-function, given respectively by 



sinh ; 



n 

k=l 



nk 



r(z)= lim ^TT(i + £ 



(Al) 
(A2) 



fc=i 



We recall the two useful formulae, holding away from 
the real negative axis, 



r(z) 

kl-"4 

r(z + x) 



2?rz 



1 

122 



+...), 



(A3) 
(A4) 



Notice the simplified asymptotic behaviors of the func- 
tions defined in Eqs. ( |l0| ) and (29) 



f(z) 



1 

127 



|x|->+oc 6Z J 







(A5) 
(A6) 



All the behavior for small arguments can be deduced 
from the one of the T-function which has a simple pole 
in the origin. In particular 
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f(z) 



1 



V2ttz 
1 

o 2z 2 



(A7) 
(A8) 



The r and functions are related to hyperbolic func- 
tions at immaginary arguments 



T(ix)T(— ix) 



2sinh7ra ' 



&(ix) + 9'(-ix) 



» sinh 7rx i 

Finally, we will need the following integrals 



C = 



dzz 4 V'(z 2 ) = 0.490 , 



dz 



( 
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0.346 



(A9) 
(A10) 

(All) 
(A12) 



APPENDIX B: LATTICE REGULARIZATION 

We present here some results about divergent quanti- 
ties computed with the lattice regularization. Using the 
lattice spectrum 



2N . fir n\ 
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and the formulae 
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it easily follows that 
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